A BLOW-UP CRITERION FOR THE 3D COMPRESSIBLE 
MAGNETOHYDRO DYNAMICS IN TERMS OF DENSITY 



ANTHONY SUEN 



Abstract. We study an initial boundary value problem for the 3D magneto- 
hydrodynamics (MHD) equations of compressible fluids in R 3 . We establish 
a blow-up criterion for the local strong solutions in terms of the density and 
magnetic field. Namely, if the density is away from vacuum (p = 0) and the 
concentration of mass (p = oo) and if the magnetic field is bounded above in 
terms of L°°-norm, then a local strong solution can be continued globally in 
time. 



1. Introduction 

We prove a blow-up criteria for the smooth solutions to the compressible mag- 
netohydrodynamics (MHD) in three space dimensions (see Cabannes [T| for a more 
comprehensive discussion on the system): 

pt + div(pu) = 0, (1.1) 

(pu j ) t + div(pu j u) + P{p) X] + (||B|% - div(B^B) = fxAu j + Xdivu Xj , (1.2) 

B\ + div(B 3 u - u ] B) = vAB j , (1.3) 

divB = 0. (1.4) 

Here u — (u^u 2 ,!* 3 ) and B = (B 1 , B 2 , B 3 ) are functions of i e i 3 and t > 
representing density, velocity and magnetic field; P — P{p) is the pressure; e, A, v 
are viscous constants. The system (|1.1[) - (|1.4[> is solved subjected to some given 
initial data: 

(p, u, B)(x, 0) = (po, u , B )(x). (1.5) 

The local existence of smooth solutions to the MHD system (|l.l[) - (ll.4j) as well 
as the global existence of smooth solutions and weak solutions are studied by many 
mathematicians in decades, see [S], [5], [T2], [TU], [2]. When the initial data is taken 
to be close to a constant state in H 3 (M. 3 ), Kawashima [TU] constructed global-in- 
time 7J 3 (K 3 )-solutions. Later, Suen and Hoff [14] generalized Kawashima's results 
to obtain global smooth solutions when the initial data is taken to be iJ 3 (R 3 ) but 
only close to a constant state in i 2 (R 3 ). The existence of global weak solutions to 
(|l.l[) - (ll.4[) with large initial data was proved by Hu and Wang [S]-[S] and Sart [T5] 
which are extensions of Lions-type weak solutions [B] for the Navier-Stokes system. 
With initial L 2 -data close to a constant state, Suen and Hoff [T3] generalized Hoff- 
type intermediate weak solutions [3]- [5] to obtain global solutions to the system 
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On the other hand, the global existence of smooth solution to the MHD system 
(| 1 . ![) - ()1 .4|) with arbitrary smooth data is still unknown. For the corresponding 
Navier-Stokes system, Z. Xin |16j proved that smooth solution will blow up in 
finite time in the whole space when the initial density has compact support, while 
Rozanova [11] showed similar results for rapidly decreasing initial density. Recently 
Fan-Jiang-Ou [2j established some blow-up criteria for the classical solutions to 3D 
compressible flows, which were further extended by Lu-Du-Yao [7] for MHD system. 

The main goal of the present paper is to generalize the corresponding results 
of Sun- Wang-Zhang [TS] to the MHD system (|1.1[1 - (|1.4[) . When the initial vacuum 
is allowed, Y. Sun, C. Wang and Z. Zhang obtained a blow-up criterion in terms 
of the upper bound of the density for the strong solution to the 3-D compressible 
Navier-Stokes equations. With the presence of magnetic field, we are able to obtain 
parallel results as in [15j except that we do not allow vacuum in the initial density. 

We now give a precise formulation of our results. First concerning the assump- 
tions on the parameters, we have 

(1.6) There exists K > such that P{p) = Kp for all p > 0; 

e, A, v > and A < e. (1.7) 

For the initial data, we assume that 

p Q ,u ,Bo £ H 3 (R 3 ) with inf(p ) > and div(S ) = (1.8) 
and we also write 

Co = \\po~ pf H s + Kilns + \\B \\ 2 Ha . (1.9) 

We make use of the following standard facts (see Ziemer [T7] Theorem 2.1.4, 
Remark 2.4.3, and Theorem 2.4.4, for example). First, given r <E [2,6] there is a 
constant C(r) such that for w £ H 1 (R 3 ), 

\\w\\ L r (R3) < C(r) (iMI^f r \\Vw\\%^] /2r ) . (1.10) 

We denote the material derivative of a given function v by v = Vt + Vi> • u, 
and if X is a Banach space we will abbreviate X 3 by X. Finally if / C [0, oo) is 
an interval, C 1 (/;X) will be the elements v £ C(I;X) such that the distribution 
derivative Vt £ X>' (K. 3 x int/) is realized as an element of C(I;X). 

We recall a local existence theorem for (jl.l[) - (|1.4[) by Kawashima [10], pg. 34-35 
and pg. 52-53: 

Theorem 1.1 (Kawashima) Assume that e, A, v are strictly positive and that the 
pressure P satisfies (1-6). Then given p > and C3 > 0, there is a positive time 
T depending on p, C3 and the parameters e, A, v, P such that if the initial data 
(po — p, uq, Bo) is given satisfying (|1.8[) and 

Co < C 3 , 

then there is a solution {p ~~u, B) to - (|1.4[) defined on M 3 x [0,T] satisfying 
p-p£ C([0,T];H 3 (M. 3 )) n ^([O, T]; iJ 2 (M 3 )) (1.11) 

and 

u,B£ C([0,T};H 3 (R 3 )) n ^([O, T]; iJ^R 3 )) n L 2 ([0, T]; iJ 4 (]R 3 )). (1.12) 
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The following is the main result of this paper: 

Theorem 1.2 Assume that the system parameters satisfy (1.6)- (|1.7p . Given p > 0, 
suppose (po — p,uq,Bq) satisfies (|1.8[) . Assume that {p — p,u,B) is the smooth 
solution as constructed in Theorem 1.1, and let T* > T be maximal existence time 
of the solution. If T* < oo, then we have 

ton, IIHU=°((0,t)xR 3 ) + HP _1 ||L = ((0,t)xR 3 ) + l|-S|U~((0,t)xR 3 ) = +°°- 

The rest of the paper is organized as follows. We begin the proofs of Theo- 
rem 1.2 in section 2 with a number of a priori bounds for local-in-time smooth 
solutions. We make an important use of estimates on the Lame operator L which 
are mainly inspired by [3] and |15) . Finally in section 3 we prove Theorem 1.2 via a 
contradiction argument by deriving higher order i? 3 -bounds for smooth solutions. 



In this section we derive a prior estimates for the local solution {p — p, u, B) on 
[0, T] with T < T* as described by Theorem 1.1. Here T* is the maximal time of 
existence which is defined in the following sense: 

Definition We call T* £ (0, oo) to be the maximal time of existence of a smooth 
solution (p-p,u,B) to (TITlT) - (fi~t|) if for any <T < T* , (p-p,u,B) solves (fi~T|) - 
(fi~4]) in [0,T]xR 3 and satisfies (TITTT]) - (fLT2"l) .- moreover, the conditions (TT^ - ^LWi 
fail to hold when T — T* . 

We will prove Theorem 1.2 using a contradiction argument. Therefore, for the 
sake of contradiction, we assume that 

l/3||i°°((0,T*)xR 3 ) + ||p _1 ||l°°((o,t*)xr 3 ) + ||-B||l=°((UT*)xR 3 ) < C. (2.1) 

To facilitate our exposition, we first define some auxiliary functionals for < 
t<T <T*: 



The following is the main theorem of this section: 

Theorem 2.1 Assume that the hypotheses and notations in Theorem 1.1 are in 
force. Given C > and p > 0, assume further that (p — p,u,B) satisfies (12.11) . 



2. A PRIOR ESTIMATES 
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Then there exists a positive number M which depends on Co, C, T* and the system 
parameters P, e, A, v such that, for < t < T < T* , 

A 1 {t)+A 2 {t)<M. (2.2) 



We prove Theorem 2.1 in a sequence of lemmas. We first derive the following 
lemma which gives estimates on the solutions of the Lame operator L = eA + (e + 
A)Vdiv. More detailed discussions can also be found in Sun- Wang-Zhang [15] . 

Lemma 2.2 Consider the following equation: 

sAv +{e + A)Vdiv(u) = J, (2.3) 

where v = (v 1 , v 2 , v 3 )(x), J = (J 1 , J 2 , J 3 )(x) with x e M 3 ande,\>0. Then for 
p € (1, oo), we have: 

(2.4) if J eW 2 -P(R 3 ), then \\D 2 x v\\ Ll , < <7||J|Up; 

(2.5) ifJ^Vcj) with (j) G W 2,P (M 3 ) ) then \\Vv\\ LP < 

(2.6) z/J = Vdiv(0) with <t> € W 2 ^(M. 3 ), then \\v\\ LP < C\\<j>\\ L P . 



Here C is a positive constant which depends only on e, X,p 

Proof. A proof can be found in [15] pg. 39 and we omit the details here. □ 
We proceed to the following a prior estimates which is the energy-balanced law: 

Lemma 2.3 Assume that the hypotheses and notations of Theorem 2.1 are in force. 
Then for any < t < T < T* , 



sup / (\p - p\ 2 + p\u\ 2 + \B\ 2 )dx + f f (|Vu| 2 + \VB\ 2 )dxds < M(C)C , 
<s<t jr3 JoJm? 

(2.4) 



where M(C) is a constant which depends on C. 
Proof. Let G = G(p) be a functional defined by 

G(p)=p [ s- 1 (P(s)~P(p))ds. 

Multiplying the momentum equation (jl.2p by u 3 , summing over j, integrating and 
making use of the continuity equation (jl.ip . we get: 



[\p\u\ 2 + G] dx 



+ f f u-div [(±\B\ 2 )I 3x3 - BB T ]dxds 



+ [ [ [e\Vu\ 2 + {e + \)(dWu) 2 ]dxds = 0. (2.5) 

JOJR 3 

Similarly, we multiply the magnetic field equation (|1.3|) by B and integrate to get 



\\B\ 2 dx 



+ 11 B ■ div (Bu T - uB T ) dxds = -v I I \VB\ 2 dxds. (2.6) 
o JoJr 3 JoJr 3 
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We then obtain (|2.4p by adding (|2.5JI to (|2.6|) and using the fact that 

/ / [u • div[(||-B| 2 )J 3X 3 - BB T }+ B ■ div(Bu T ~uB T )]dxds = 0. 

JoJr 3 

We obtain the following L A bounds for u and B: 



□ 



Lemma 2.4 Assume that the hypotheses and notations o/Theorem 2.1 are in force. 
Then for any < t < T < T* , 

(\u{x,t)\ 4 + \B(x,t)\ 4 )dx < M. (2.7) 

j 

Proof. Multiply (|1.2I) by 2|w| 2 w and integrate to obtain 

/ p|u| 4 da: + / 2|ii| 2 [e|VM| 2 + (A + e)(div(u)) 2 ]da; 
Jm. 3 Jr 3 

+ [ 8 [e\u\ 2 \V (\u\)\ 2 + (e + X)(div (u))\u\u -V{\u\)] dx 
Jr 3 

{P{p) - P{p))d\v{\u\ 2 u)dx 

+ j 2\B\ 2 div(\u\ 2 u)dx + / A\u\ 2 u ■ div{BB T )dx (2.8) 
Jr 3 Jr 3 



The third term on the left side of (|2.8[) can be estimated from below by 
8[e|u| 2 |V(|u|)| 2 + (e + A)(div(«))|u|« • V(|u|)]da; 



> 



4\u\ 



e|Vu| 2 + 2(e- £±^)|V(M)| 2 



dx. 



By assumption (|1.7p we have e < A, hence it implies 



\u\ 2 (s 



e + A 



)|V(|m|)| 2 > M \ \u\ 2 \Vu\ 2 dx. 



(2.9) 



On the other hand, we multiply (jl.3j) by 4\B\ 2 B and integrate to get 
/ |S| 2 da;+ / 4is\B\ 2 \\7B\ 2 dx+ [ 2v\B\ 2 \\7B\ 2 dx 

dt JR 3 JR3 J R 3 

= - I \B\ 2 B ■ div(Bu T - uB T )dx. 
Jr 3 

Adding (12.10)) to (|2.9I) and integrate with respect to t, we get 

(p\u\ 4 + \B\ 4 )dx] + [ [ (\u\ 2 \Vu\ 2 + \B\ 2 \VB\ 2 )dxds 
J JoJr 3 

i(P(p) - P(p))div(\u\ 2 u)dxds + [ [ 2\B\ 2 div(\u\ 2 u)dxds 

JoJr 3 

[ [ 2\u\ 2 u-div(BB T )dxds + [ [ \B\ 2 B • div{Bu T - uB 1 
JoJr 3 JoJr 3 



(2.10) 



< M 



"I 

lOJR 3 



)dxds . 
(2.11) 
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Using the assumption (|2.1|) . the right side of (|2.1ip can be bounded by 



(\P(p)-P(p)\ 2 \u\ 2 + \B\*\u\ 2 )dxd S 



< M 



T*C 



\B\ 



ojr 3 



\u\ )dxds 



I I (|u| 2 |Vu| 2 + \B\ 2 \VB\ 2 )dxds 
JoJr 3 



(|u| 2 |Vw| 2 + \B\ 2 \VB\ 2 )dxds 



OJR 3 



(2.12) 



Using (I2.12p on (|2.11l) and applying Cauchy Inequality, we get 

/ (\u\ 4 + \B\ 4 )dx < M + // (|u| 4 + \B\ 4 )dxds, 
Jr 3 JqJr 3 

and (|2.7p now follows by Gronwall's inequality. □ 
We obtain estimates on the functional A\ in terms of H: 

Lemma 2.5 Assume that the hypotheses and notations o/Theorem 2.1 are in force. 
Then for any < t < T < T* , 

Ai(t) < M[l + H(t)]. (2.13) 
Proof. We multiply (ll.2[) by v? , sum over j and integrate to get 

|Vm| 2 c&e+ / / p\u\ 2 dxds 
JoJr 3 



[u.V(hB\ 2 )-u-div(BB T )] 

OJR 3 L 



OJR 3 



\Vu\ 3 dxds. (2.14) 



Niext we multiply (| 1 .3[) by B t and integrate, 



\VB\ 2 dx 



\B t \ 2 dxds < C 



OJR 3 



B t ■ div(uB T - u T B)dxds 



(2.15) 



Adding (|2~14|) and (|2~T5]) . we obtain 



/ (|Vu| 2 + |VB| 2 )dx+ / / (p\u\ 2 + \B t \ 2 )dxds 

Jr 3 JoJr 3 



<C + 



[ [ \Vu\ 3 dxds+ f f {\VB\ 2 \B\ 2 + \Vu\ 2 \B\ 2 + \VB\ 2 \u\ 2 )dxds. 
JoJr 3 JoJr 3 



(2.16) 



The second term on the right side of (|2.16p is bounded by 



\\7u\ 2 dxds 



OJR 3 



OJR 3 



\Vu\ 4 dxds) <C Q + H(t), 



where the last inequality follows by Lemma 2.3. For the last integral on the right 
side of (|2.16p . using assumption (|2.1[) . it can be bounded by / / \VB\ 2 \u\ 2 dxds + 



M 



t JOJR 3 

/ / (|Vu| 2 + \VB\ 2 )dxds. So it remains to estimate / / V B\ 2 \u\ 2 dxds. 
J oJr 3 JoJr 3 



Recall from Lemma 2.4 that, for < t < T < T* 

\u\ 4 dx < M, 
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Therefore, using (|1.10j) . 



VB\ 2 \u\ 2 dxds < 



qjr 3 



< M 



< M 



< MCI 



< MCI 



\7B\ 4 dx 
VB\ 2 dx 



\VB\ 2 dxds 



dx ds 



\DiB\ 2 dx ds 



QJR 3 



\AB\ 2 dxds 



qjr 3 



qjr 3 



{\B t \ 2 + \VB\ 2 \u\ 2 + \Vu\ 2 \B\ 2 )dxds 



A 1 {t) + C Q + 



VB\ 2 \u\ 2 dxds 



and by Cauchy inequality, we obtain 



QJR 3 



VB\ 2 \u\ 2 dxds < M 1 + 4(4)' 



Applying (pTrTj) to (EHo) and absorbing terms, (|2TT3)) follows. 



(2.17) 
□ 



We derive the following estimates on the effective viscous flux which were first 
described by Hoff [3] and later modified by Sun- Wang-Zhang [15] . 

Lemma 2.6 Assume that the hypotheses and notations o/Theorem 2.1 are in force. 
Then for any < t < T < T* , 



sup 



IVwI 



w t \ 2 dxds+ / / \D 2 w\ 2 dxds < M, (2.18) 
JqJr 3 



0<B<tJM3 

where w = u — v with v satisfying: 

eAv + (e + A)Vdiv(v) = V(P(p) - P{p)). 

Proof. Using the momentum eqaution (jl.2l) . 

pw t - sAw - (e + A)Vdiv(w) = -pu ■ Vu - pv t - V(i|5| 2 ) + div(BB T ). (2.19) 

Multiply (12.19)) by w t and integrate, 



elVwl dx 



o Jo 



(e + A) | divw \ 2 dxds - 



OJR 3 



p\wt\ 2 dxds 



[-pu ■ Vm - pv t - V(-\B\ 2 ) + diviBB 1 )] • w t dxds. (2.20) 
qJr 3 2 

The first term on the right side of (|2.20[) can be estimated as follows: 



OJR 3 



(— pu ■ Vit) • w t dxds < M 



\u\ 2 \Vu\ 2 dxds 



Uo 



\w t \ 2 dxds 



QJE. 3 



(2.21) 
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For the term dxds, using Lemma 2.4, 

JoJR 3 



\u\ 2 \\Ju\ 2 dxds < 



QJR 3 



\u\ dx 



o VE 3 



|Vu| dx 



< M 



< M 



|Vw| 4 + \\7v\ 4 )dx ) ds 



VR 3 



\Vw\ 2 dx 



\D„w\ dx I ds 



M 



o \Jm 3 



\Vv\ 4 dx 



< 



M ( f [ \Vw\ 2 dxds) ( [ [ \D 2 x w\ 2 dxds] 
\JoJr 3 J \Jojr 3 J 

\P(p)-P(p)\ 4 dx) ds 



< M 



M 



Wwrdxds 



WR 3 



\D rr .w\ dxds 



MC£T*K 



(2.22) 

where the last inequality follows by Lemma 2.2 and assumption (I2.1[) . Therefore 
(|2T2T|) becomes 
ft f 

{—pu ■ Vu) • Wtdxds 



< M 



\Vw\ dxds 



ojr 3 



-l2„..|2. 





i 

2 
















dxds 


Jqjr 3 





(2.23) 



The third and the fourth term on the right side of (|2.20p are bounded by 



\V B\\B\\w t \dxds < M 



< MCI 



\VB\ 2 \B\ 2 dxds 



J R ;1 



\wA 2 dxds 



\wA dxds 



ojr 3 



(2.24) 



It remains to estimate the term / / — pvt ■ Wtdxds on the right side of (|2.20p . By 

JoJr 3 

the definition of v and P{p), we have 

eAv t + (e + A)Vdiv(u t ) = VP{p) t = Vdiv(-P(p)u), 
Hence we can apply Lemma 2.2 and Lemma 2.3 to get 

v t \ 2 dxds< [ [ \P(p)u\ 2 dxds < MC a , 



OJR 3 
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and 



-put ■ Wtdxds < M 



U! t 



| dxds I 



\v t \ 2 dxds 



<M[ / / \w t \ 2 dxds 



OJR 3 



(2.25) 



Using (|2~2"5j) . (j2~M| and (g^ on (|2~2Ut , 



e|Vw;| 2 da; 



->o 



+ / / (e + A)|div(w)| 2 efeds + / / p\w t \ dxds 



ojm. 3 



< M { \Vw\ 2 dxds 



4 / ft 



D 2 x w\ 2 dxds) +M. (2.26) 



0./R 3 



It remains to estimate / / \D 2 w\ 2 dxds. We rearrange the terms in (|2.19|) to get 



eAw + (e + A)Vdiv(w) = pw 4 + pVit • u + pv t + V(-|5| 2 ) - div(BB T ) 
and so by Lemma 2.2, 



/o 

< M 

Therefore 



ojr 3 



D 2 w\ 2 dxds< / / (\pw t \ 2 + \pVu-u\ 2 + \pv t \ 2 + \VB\ 2 \B\ 2 )dxds 



\Dlw\ 2 dxds + 1 



Wt\ 2 dxds + ( / / Vwl cfods 



/o 



\wt\ 2 dxds 



OJR 3 



\Vw\ 2 dxds + 1 



and we apply the above to (|2.26p to conclude 



etywrdx 



+ 11 (e + X)\div{w)\ 2 dxds + / / p\w t \ 2 dxds 



JoJR 3 









< M 


[fl 


\Vw\ 2 dxds + 1 




J OJR 


3 



which (|2.18|) follows by Gronwall's inequality. 



□ 



We finally obtain an estimate on the functional Ai which is sufficient to prove 
Theorem 2.1: 

Lemma 2.7 Assume that the hypotheses and notations o/Theorem 2.1 are in force. 
Then for any < t < T < T* , 



A 2 (t)<M[A 1 (t)+H(t) + l] 



(2.27) 
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Proof. Taking the convective derivative in the momentum equation (|1.2[) , multiply- 
ing it by u J , summing over j and integrating, 

rt 



2 dx + (j \Vii\ 2 dxds 
JoJr 3 



<C + H(t)+ \B\ 2 (\B t \ 2 + \u\ 2 \VB\ 2 )dxds. (2.28) 



sup 

0<s<t 



Next we differentiate the magnetic field equation (|1.3[) with respect to t, multiply 
by B t and integrate, 

t 



\B t \ 2 dx 



' I I \VB t \ 2 dxds 

JoJr 3 



B t ■ [div(Bu T - uB T )] t dxds. 
Adding the above to (|2.28[) and absorbing terms, 



sup 

0<s<t 



\B t \ 2 )dx+ I I (|Vu| 2 + \VB t \ 2 )dxds 
JoJr 3 



< M 



Ai+ H 



M 



ojr 3 



\B\ 2 \u\ z {\Vu\ z + \VB\ 2 )dxds 



ojr 3 



(\B\ 2 \B t \ 2 + \B\ 2 \u\ 2 + \B t \ 2 \u\ 2 )dxds (2.29) 



The third term on the right side of (|2.29p is bounded by 



|tt| 2 (|Vu| 2 + \VB\ 2 )dxds < 



\u\ |Vu| dxds 



< M 



OJR 3 

l + Aj 



OJR 3 



\u\ 2 \VB\ 2 dxds 



(2.30) 



where the last inequality follows by (I2.17|) and (|2.22l) . The last term on the right 
side of (|2.29p is bounded by 



/ / (\Bf + \ii\ 2 )dxds + I I \B t \ 2 \u\ 2 dxds 
JoJr 3 JoJr 3 



< MAi 



\JR 3 
t 



< MAi + M 



< MA! + M 



\JR 3 



IB^dx 



\B t \ 2 dx 



\u\ dx I ds 



\VB t \ 2 dx) ds 



\B*\ 2 dxds 



QJR 3 



OJR 3 



\VB t \ 2 dxds 



< MA\ + A{ A* . 

Using ([2~3D|i and (|2~3T1) on ([2~2^|l and absorbing terms, ([2~2"g| follows. 

proof of Theorem 2.1. Recall from Lemma 2.5 and 2.7 that 

A x < M [H + 1] 



(2.31) 
□ 
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and 



A 2 <M[H + A x ] 



So it remains to estimate H. Let w and v be as denned in Lemma 2.6. Then 



Vitl dxds < 



I Viol dxds 



OJR 3 



|Vw| dxds 



OJR. 3 



(2.32) 



The second term on the right side of (|2.32l) is bounded by / / \P(p) — P(p)\ 4 dxds. 
And for [ [ \Vw\ 4 dxds, using ([235]) . 



|Viy| <ia;<is < 



VR 3 



< ( sup / \D%w\ 2 dx I I sup / \S7w\ 2 dx 



0<s<t 



< M [ sup / |D 2 w| 2 da; 

V0<s<t 



0<s<t 



Notice that, by rearranging the terms in (|2.19l) . 

eAw + (e + A)Vdiv(w) = pii + V(i|B| 2 ) - div{BB T ), 

and so by Lemma 2.2, 

/ \D 2 x w\ 2 dx < [ \pu\ 2 dx+ [ \V{l-\B\ 2 )~div{BB T )\ 2 dx 
Js. 3 Jr 3 Jr 3 2 

<M I \ii\ 2 dx+ [ \\7B\ 2 \B\ 2 dx 
Jr 3 Jr 3 

< M[A 2 + At]. 
Therefore we conclude that 

H <M[A 2 +Ai]? , 

and (I2~2l follows. 



\Vw\ 2 dx 



OJR 3 



□ 



3. Higher Order Estimates and proof of Theorem 1.2 



In this section we continue to obtain higher order estimates on the smooth local 
solution (p — p, u, B) as described in section 2. Together with Theorem 2.1, we show 
that, under the assumption (|2.1[) . the smooth local solution to (ll.l[) - (|1.4l) can be 
extended beyond the maximal time of existence T* as defined in section 2, thereby 
contradicting the maximality of T*. The following is the main theorem of this 
section: 

Theorem 3.1 Assume that the hypotheses and notations in Theorem 2.1 are in 
force. Given C > and p > 0, assume further that (p — p,u,B) satisfies (|2.1| . 
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Then there exists a positive number M' which depends on Co, C,T* and the system 
parameters P, e, A, v such that, for < t < T < T* , 

sup \\{p-p,u,B)\\ H s m + [ \\(u, B)(; s)\\ 2 HHRa) ds <M' (3.1) 

0<s<t Jo 

Proof. We give the proof in a sequence of steps. Most of the details are reminiscent 
of Suen and Hoff [14j and we omit those which are identical to or nearly identical to 
arguments given in [14] . We first begin with the following estimates on the effective 
viscous flux F and the vorticity matrix uj: 

Step 1: Define 

F=(2s + X)dW(u)-(P(p)-P(p)), 
u = u>i> k =u{ k -u k x ., 

Then for q € (1, oo), 

||Vtt(-,i)||x. < M(q) [\\F(;t)\\ L « + \M;t)\\L« + \\(P( P ) P(p))(;t)\\^} , (3.2) 

||Vw(- ) t)||i«<(«)M[||pu(- ) t)|| i « + ||VB.B(-,t)|| £ ,], (3.3) 
where M(q) is a positive constant depending on q and 



sup / (\VF\ 2 + \\7oj\ 2 )dx < M'. 

0<s<t JR3 



(3.4) 



proof of Step 1. We give the proof of (|3.2[) as an example. Using the definition of 
F and to, 

(2e + X)Au j = F Xj + (2e + A)< fe + (P(p) - P{~p) X] , (3.5) 
Differentiating and taking the Fourier transform we then obtain 

(2s + \)u x (y,t) = f^F(y, t) + (2s + X)f^^(y, t) + f^(f^P)(y, t) 

\v\ \v\ \y\ 

and (I3.2p then follows immediately from the Marcinkiewicz multiplier theorem 
(Stein 13; , pg. 96). Similarly, (I3.3[) can be proved by the same method. Also, 
by the definition of F, we have 

AF = div( 3 ), (3.6) 
where — pu? + (\\B\ 2 ) Xj — div(S : 'i?). So we have 



sup / \VF\ 2 dx < sup / \g\ 2 dx 

0<s<tJR 3 0<s<tjR3 

< sup f (p\ii\ 2 + \VB\ 2 \B\ 2 )dx < M', 

0<s<t JR3 

and similarly, sup / \Vtu\ 2 dx < M' , which proves (|3.4p . □ 

0<s<t JR3 

Step 2: The velocity gradient satisfies the following bound 

\\\7u(-,t)\\ L °ads < M'. 



proof of Step 2. The proof is identical to Suen and Hoff [2] pg. 51-53, and we 
omit the details here. □ 
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Step 3: We further obtain 

||2^u(-,t)|| L 2 < M' \\\fm(;t)\\L> + IIVB • B(-,t)\\ L 2 + ||VP(-,t)|| L2 



(3.7) 



\\D 3 x u(;t)\\ L 2 < M' [||Vp-ti(.,i)|| za + ||pV«(-,t)|| La + H-B- D 2 x B(;t)\\ L2 ] 

+ M' [|||VS| 2 (.,i)|| ia + \\D 2 x P(.,t)\\ L *] . (3.8) 



proof of Step 3. These follow immediately from the momentum equation (|1.2j) and 
the ellipticity of the Lame operator eA + (e + A)Vdiv. □ 

Step 4: The following H 2 -bound for density holds 



sup \\(p~ p)(-,s)\\ H 2 < M' . 

0<s<t 



(3.9) 



proof of Step 4- We take the spatial gradient of the mass equation multiply 
by V/9 and integrate by parts to obtain 



From (RTF] 



^ / \Vp\ 2 dx<M' f \Vp\ 2 dx+ [ \D 2 x u\ 2 dx 

JR3 [Jr3 J R 3 



D x u\ 2 dxds< / / (\u\ 2 + \VB ■ B\ 2 + \\7p\ 2 )dxds 



(3.10) 



QJR 3 



< 



M'+ [ [ \Vp\ 2 dxds. 



Applying the above to (|3.10|) and using the result of Step 2, 

sup \\Vp(-,s)\\ L 2 < M 1 . 

0<s<t 

By similar argument, we can show that sup \\D 2 p(-, s)||^2 < M' and (|3.9|) follows. 

0<s<t 

□ 

Step 5: The velocity and magnetic field satisfy 

sup (||u(-,«)|| fl » + ||S(- ) fl)|| H 8)<M / . (3.11) 

0<s<t 

proof of Step 5. Define the forward difference of quotient D' t l by 

Dl(f)(t) = (f(t + h)-f(t))h- 1 

and let = D^(yP) + u ■ Vu J . By differentiating the momentum equation, we 
obtain 

/ p\E X] \ 2 dx+ f f (\VE Xj \ 2 + \D?(div(u x .) + u- V(div(u X] )\ 2 )dxds 

JR 3 JoJR 3 

M'+ [ [ \\7E\ 2 dxds + 0(h), 

JCWR 3 

where 0(h) — > as h — > 0. Therefore by taking h — > 0, 

sup ||Vu(-,s)|| L 2 + f j \D 2 x u\ 2 dxds < M' . 

0<s<t JoJR 3 

The bound for V B t can be derived in an exactly same way. □ 



< 
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Step 6: Finally we have the following bounds 



[ [ (\D 4 x u\ 2 + \D 4 x B\ 2 )dxds < M' 1+ f f \D x p\ 2 dxds 
JoJr 3 JoJtsl 3 



(3.12) 



sup {\\D 3 x p(; S )\\ L 2 + \\DlB(; S )\\ L 2) + I I \D x u\ 2 dxds<M'. (3.13) 

0<s<t 



proof of Step 6. For (|3. 12|) , it can be obtained by differentiating (11.2[) and (II .3[) 
twice with respect to space, expressing the fourth derivatives of u and B in the 
terms second derivatives of u, Bt, \7p and lower order terms, and applying the 
bounds in ([2~Tj) and (|3~TT]) 

For (|3.13[) . it can be obtained by applying two space derivatives and one spatial 
difference operator D x . defined by 

D h Xj (f)(t) = (f(x + he J )-f(x))h- 1 

such that 

/ \D h Xj D Xi D XkP \ 2 dx < W + I I (\D x u\ 2 + \D h Xj D Xi D XkP \ 2 )dxds 
Jr 3 JoJr 3 

<M'+ [ [ \D 3 x p\ 2 dxds. 

Taking h — > and applying Gronwall's inequality, we obtain the required bound 
for the term \\D x p(-, s)\\ L 2, □ 

proof of Theorem 1.2. Using Theorem 3.1, we can apply an open-closed argument 
on the time interval which is identical to the one given in Suen and Hoff |14j 
pp. 31 to extend the local solution (p — p,u,B) beyond T* , which contradicts the 
maximality of T* . Therefore the assumption (|2.ip does not hold and this completes 
the proof of Theorem 1.2. □ 
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